FIG. 1 
FIG. 2 

Al: Downloading intake air amount Qa and fuel injection amount Qb 

al: Calculating excess air ratio X , and amount Me of 

particulates discharged in accordance with excess air ratio X 

X vs. Me in map ml 

A2: Calculating Mb of particulates burnt per unit time t 
(Mb=a XPMXt) 

A3: 

Calculating total Mampt of Ma of particulates accumulated 
per unit time based on Me, Mb & Ma of discharged, burnt and 
accumulated particulates 

Mampt > predetermined value? 

Executing forced regeneration 

Catalyst temperature gt (=filter temperattire) 
bO: Calculating fuel biurning velocity coefficient a 

a vs. catalyst temperature map mO 
Fig. 3(a) Large Small Large 
Fig. 3(b) Large Small Small 
Fig. 3(c) Large Small Low High 
Fig, 4(a) 
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Determination of excess air ratio frequency Time 



Fig. 4(b) 

Excess air ratio frequency y 
High Low 

Fig. 5(a) 

Torque (fuel injection amount) 
Large Large Small 
Engine speed High 

Fig. 5(b) Large 

Fig. 6 

Forced regeneration routine 
SI: Calculating Me 
S2: Calculating Mb 
S3: Calculating Ma 
S4: Ma > Ma a 

S5: Executing forced regeneration 
Return 

Fig, 7 Crank angle 
Fig. 8 

Downloading intake air amount Qa and fuel injection amount Qb 

al': Calctdating excess air ratio X 

a2-l': 

Calculating excess air ratio frequency y A t, 

and setting 1 when excess air ratio X is equal to or 

less than predetermined value (otherwise setting 0) 

a2-2': 

Totaling Ma At : Ma At =f( y At) 
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bl: 

Calculating filter temperature frequency i3 At, and setting 1 when 
temperature is equal to or less than the predetermined value (otherwise 
setting 0) 

NOx/Soot ratio 

A2': 

Calculating particulate burning coefficient a A t (= f( jS A t) , 

A3": 

Calculating Mei using 

PMi = PMi-i + (MaA MbAt)XAt 

Ma ^ predetermined value? 

Executing forced regeneration 

Catalyst temperature gt (=filter temperature) 

b4": 

Calculating Mb A t: (Mb A t = a A t X PMi-i) 
Fig. 9(a) 

Forced regeneration (timing detection) routine 
SIO: Calculating Ma A t 
S20: Calculating Mb A t 

S30: Calculating PMi (= PMi-i + (Ma At-MbAt)*At) 
S40: PMi ^ predetermined value? 
S50: Executing forced regeneration 
Return 

Fig. 9(b) 

Calculating M a At 



25 



Sll: Downloading Qa and Qf 

S12: Calculating excess air ratio X 

S13: Calculating excess air ratio frequency y At 

S14: Calculating Ma A t(=f ( v A t) 

End 

Fig, 9(c) 

Calculating Mb A t 

s21: Downloading catalyst temperate gt 

s22: Calculating filter temperature frequency jS A t 

s23: Calculating particulate burning velocity coefficient a A t 

s24: Calculating MbAt(=aAtX PMi-i) 

End 



